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Abstract 

The flavour dependence of the mixed quark-gluon condensate is studied through 
the analysis of correlators of the hybrid current a p = gs'y P ^G Pll d. The flavour 
symmetry breaking for this type of condensates is found to be less than that for 
the quark condensates. For the ratio of strange to nonstrange condensates we 
obtain R = 0.95 ± 0.15. For the kaon coupling to the current a p we find 5' 2 = 
(0.020 ± 0.005)Gey 2 , which is an order of magnitude smaller than analogous chi- 
rally unsuppressed coupling. 
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Introduction 



The mixed quark gluon condensate (0|scr M1/ (jr M! ,s|0) is an important phenomenological 
parameter characterizing the structure of QCD vacuum. Considerable attention has been 
paid in the past to determination of it's value, as well as of the ratio 

(da fj, v G fii/d) 

of strange and nonstrange quark-gluon condensates [1-5]. 

A QCD sum rule (SR) analysis of the baryon decuplet indicates the value R ~ 1.3 || 
while various hybrid sum rules favour the values in the lower range R — 0.5 — 0.85 1], 0- • 
At the same time, the analysis of heavy to light mesons || results in the restriction 
1 -R< 0.15. 

The two point functions of the hybrid current a M = gsjpj^Gp^d which have been 
studied in (I], |2[ by the QCD sum rules technique are potentially the most sensitive to 
the value of R. The reason is, that in this case the leading contributions come from the 
mixed quark-gluon condensates of interest and the correlators are themselves proportional 
to the SU(3) breaking. In the present work we continue the study of the correlators of 
the hybrid current a M , being motivated by several reasons: 

First, we are interested in the value of the kaon coupling to the hybrid current 

(0\a„\K(p)) = -SkPpp 

which was left out in the the previous analyses. This parameter measures the Chiral 
symmetry breaking for kaon (in the chiral limit i5' 2 = 0). Besides, this quantity can 
contribute significantly to various kaon transitions calculated in the Vertex Sum Rules 
approach (see e.g. [0]). 

Next, we have calculated some additional contributions to theoretical expressions 
which happen to be important numerically and affect the value of the strange quark- 
gluon condensate. The parameter set which we use also differs from that of previous 
estimates. 

For the determination of the unknown qantities we apply the recently proposed proce- 
dure of simultaneous analysis of several sum rules || which we believe provides the means 
of more detailed study for the type of problems considered here. The allowed range of un- 
known parameters is determined through the requirement of the compatibility of different 
SRs which is tested using the \ 2 criterium. 

Equations 

Consider the correlation functions of a M with pseudoscalar and axial vector currents: 

ig J dxe^(0\T (d W s(x), s W G pu d(0)) |0> = g^TS? (q 2 ) + q^q v Uf (q 2 ) (2) 

ig J dxe^ x (0\T (d lB s(x), s lpl5 G ptl d(0)) |0> = q,W 3 {q 2 ) (3) 
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The operator product expansion for the invariant functions in (2), (3) is obtained in 
the Euclidean region by calculating diagrams shown in fig.l and has the following form: 



n-(-Q 2 ) = g^Lzfegg) + j£_ g2 (in 2 ( Q 2 /p 2 ) - HQ 2 /, 2 )) 

«3d> - (ss)) HQ 2 /, 2 ) + ms{a ^f 2) {HQ 2 /, 2 ) ~ l) 
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Here we work in the chiral limit = for the d quark and restrict ourselves to the 
terms linear in SU(3) breaking. The calculation has been performed in the MS scheme 
and the fixed point gauge technique has been used for condensate contributions. 

The quark and gluon condensate contributions are given by the one loop diagrams 
b),c) and d). The latter one requires the removal of infrared divergence which occurs due 
to the mixing of the operators m s G 2 and igsGs. In spite of their one loop suppression, the 
quark and gluon condensate contributions are quite important and comparable to that of 
the quark-gluon condensates. Inclusion of these terms shifts R to its higher values, but 
mostly affects 5' 2 in comparison with ||. 

In (4-6) the four quark condensates are reduced to the squares of quark condensates, 
but we retain a factor p to account for the violation of vacuum dominance. The tree level 
dimension seven contributions have been evaluated in the vacuum saturation approxima- 
tion. However, due to numerical smallness of these terms the approximation has only 
minor effect on our final results. 

We keep the perturbative contribution (diagram a) ) only for U ps where it is of the 
order 0(m s ) , for the similar contribution is proportional to m 2 and thus negligible 
in our approximation. The gluon condensate does not contribute to equations (5), (6) to 
the considered order, since it is also 0{m 2 ) p|. 

The anomalous dimensions of the operators are accounted for by taking e.g. 



( f . G9 )(,)H^)(l G en(^P) 3,/54 ; m ,^l^(^p)~ 4/S (7) 

We saturate the phenomenological sides of SRs by the lowest lying intermediate states: 
the K meson for the pseudoscalar correlator, and the K and K\ for the axial vector one. 
We use the usual model spectra with K and K\ as narrow resonances and the continuum, 
equal to the theoretical one, starting at some threshold s . Following 0, instead of 
the two nearby resonances i^i(1270) and i^i(1400) we substitute an effective resonance 
Ka(1335). 

After the Borel transformation || we obtain the following set of equations: 

^ l - R )^Gd) + a, , _ . _ e _ XQ) 



AM 2 3tt 
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where C is the residue of the K\ effective pole, and xq = jfe with s being the continuum 
onset. The running parameters are taken at fi 2 = M 2 . 

Eqs.(8-10) constitute a system of sum rules to be used for the determination of the un- 
known quantities R, 5' 2 and C. Following the previous analyses |2| we could eliminate 
the couplings C and 5' 2 from these equations and express R as a function of QCD param- 
eters and M 2 . This however would mean that the two of three equations are required to 
hold exactly at any M 2 , while the third one is considered as approximate equation with 
all uncertainties and errors accumulated in it. Instead, we follow the procedure developed 
in || and consider the Borel transformed SRs (8-10) as a system of independent approx- 
imate equations, where we try to estimate the expected errors and find a set of unknown 
parameters which makes the whole system maximally consistent. 

The error analysis of such a system is the most problematic task where the model 
assumptions are unavoidable. Following || we model the errors of equations (8-10) in 
the following way: For each one of these equations we choose some reference point M 2 
(i=8,9,10) of the Borel parameter, where the relative error Wi is assumed to be minimal. 
Then the absolute error distribution is described by a simple function 

D(M 2 )J^ *M 2 <M 2 GeV 2 
[1 HM 2 >M 2 GeV 2 

which has a pole at Mq < M 2 (to account for the divergence of the power series) and is 
constant for M 2 > M 2 . 

Finally the absolute error is calculated as 

\(M 2 ) =D{M 2 )u l A l (12) 

where A$ = Ai(M 2 ) denote the r.h.s.'s of eqs.(8-10) at the reference values of the Borel 
parameter. Note, that due to the decrease of Ai(M 2 ) at high M 2 eq.fll2"[) corresponds to 
growth of the relative errors for both high and low values of this parameter. 

To get an idea of the scale of expected errors, we first examine each of the equations 
(8-10) separately and check their stability with respect to M 2 for different values of (sGs). 
Fig. 2 shows the M 2 dependence of d~' 2 and C (normalized to their final values) as defined 
from eqs.(8-10) with R = 0.85. The SR (El) is the most sensitive to the value of R. It 
becomes more stable for higher values of this parameter, but as can be seen from the 
figure, even in this case the stability is not very good, which can be attributed to the 
unknown higher orders corrections or the roughness of the model of the phenomenological 
spectrum. On the other hand eq.@ is the most stable in M 2 for a wide range of R and 
we consider it to be the most reliable one. For eq.fllOD the stability plateau starts at 
higher M 2 values than for eqs.(9-10), which can be considered as indication of the higher 
characteristic mass scale for this equation [yr? K ~ 1.8GeV 2 ). Thus our check shows that 
we should expect a greater relative error in (pj) than in other equations and use higher 
M 2 values for SR (10). 
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Results and Discussion 



In order to fix the unknown quantities R, 8 and C we first scan all their possible values. 
For each set we evaluate the equations (8-10) at four separated M 2 points, within the range 
of their expected validity. Thus we build up a system of 12 linear approximate equations 
where for each one we calculate the errors according to eqs. (11), (12). We look for a set 
of unknown parameters which makes these 12 equations maximally consistent, e.i. gives 
minimal Xdof- Once we find such a solution, we test its stability against variations of 
different parameters involved in the derivation (M 2 points, Ui etc.) and adjust the latter 
to get the least sensitivite result. Note, that this is not a standard x 2 procedure since in 
our case the errors depend on the solutions themselves. The details and justification of 
the approach can be found in ||, where we have used the similar method to restrict the 
values of standard condensates. 

For eqs. (8), (9) we take the M 2 points to be in the range (0.8 4- 2) GeV 2 , the effective 
pole in the error distribution is set to M = 0.7GeV 2 and we use M 2 = 1.2GeV 2 as a 
reference value of M 2 . For eq.(ll) we increase all these values by (0.2 — 0A)GeV 2 . The 
typical relative errors which we have used in our analysis are: cjg — 0.25, u> 9 = 0.15 and 
uj 10 = 0.20. This is quite a conservative choice, especially taking into account that the 
relative errors grow for both higher and lower values of M 2 . The continuum onset in 
equations (9), (10) is taken as s ~ 2.5GeV 2 , while in the equation (||) where there is no 
K\ contribution, the continuum threshold is taken to be lower, so — 1.5GeV 2 . 

For the numeric evaluation we use the following values of the input parameters: 

Ams = 150MeV, m s = 180MeV[ll, 0, 
(qq) = -(0.24Gev) 3 , (ss) = (0.7 ± 0.1) (13) 

(-G 2 ) = 1.2 • 10- 2 GeV 4 §, (dGd)/(dd) = OMGeV 2 , p = 4@.F] 

7T 

The minimal x 2 solution is given by R = 0.93, 5' 2 = 0.021GeV 2 and C = -6 ■ 10~ 4 GeV 4 
corresponding to Xmin = 0-6- This solution is quite stable with respect to the varia- 
tions of different parameters introduced in our procedure. The changes of M 2 points 
by ±0.2GeV 2 , of w^s by 20% and M by ±0.1GeV 2 do not affect the result signifi- 
cantly. We have also checked the dependence on the parameters of eq.(13). So, changing 
(ss)/(dd) = 0.6(0.8) we get R = 0.94(0.89), 5' 2 = 0.02(0.023); m s = 160(200)Me1/ gives 
R = 0.92(0.95), 5' 2 = 0.019(0.024)GW 2 and M 2 = 0.8(0.5)Ge1/ 2 gives R = 0.91(0.96) 
and 5' 2 = 0. 024(0. 019)GeV 2 respectively. The final result is also insensitive to the varia- 
tions of s and s' . 

In fig. 3 we have plotted the ranges in the (R, <5 /2 ) plane corresponding to the solutions 
with Xd.o.f < 1 an d Xd.oj < Xmin + 1 = 1-6- We consider all the solutions having 
Xd.o.f < Xmin + 1 = 1-6 as acceptable and thus quote our final result as 

# = 0.95 ±0.15 5' 2 = (0.020 ±0.005)GeV 2 (14) 

From the figure one can notice that the correlation of these quantities is insignificant. 

lr This corresponds to the results of Qwhere instead of the present value we have used {y/alqq) = 
-(0.24Gev) 3 of §. 
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The value of R is larger than that in previous estimates [1-3] and is consistent with 
the symmetry limit R — 1. The errors are also bigger, but to our feeling here they look 
more realistic. The result conforms well with the conjecture of ||, based on the analysis 
of heavy-to-light systems, that inclusion of the gluon field into the operator reduces the 
strength of the flavor symmetry violation. 

The parameter 5' 2 is an order of magnitude smaller than the similar, chirally unsup- 
pressed coupling of the kaon to another hybrid current s^ v G Vil d (5 2 ~ 0.2GeV 2 
Consequently, it will not contribute significantly to the vertex sum rules except the cases 
proportional to the SU(3) flavour symmetry breaking. 
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Figure captions 



Fig.2 Stability test for SRs (8-10) with R = 0.85: 
• — 5 >2 /(0.002GeV 2 ) from eq.(8); 

x — <5 /2 /(0.002G'e\/ 2 ) from eq.(10) with C = 6 • 10- 4 GeV 4 
+ - C/(-6 • 10- 4 GeV 4 ) determined from eq.(9). 



Fig. 3 The allowed region in the plane (R, S' 2 ). The areas with Xd.o.f < 1 an d 
xlo.f < Xmin + 1 are indicated. 
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